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Abstract 

The Bethe ansatz equations of the fundamental Sp(2N) integrable model are solved 
by a peculiar configuration of roots leading us to determine the nature of the excitations. 
They consist of iV elementary generalized spinons and N — 1 composite excitations made 
by special convolutions between the spinons. This fact is essential to determine the low- 
energy behaviour which is argued to be described in terms of 2N Majorana fermions. 
Our results have practical applications to spin-orbital systems and also shed new light to 
the connection between integrable models and Wess-Zumino-Witten field theories. 
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The study of quantum one-dimensional integrable models has turned out to be a fruitful 
venture since the seminal work of Bethe in 1931 [|I|]. Over the years, solvable models have been 
extremely useful in many subfields of physics, providing us a rich laboratory in which new 
theoretical insights and non-perturbative methods can readily be tested. Recent progress in 
the experimental study of low- dimensional materials, e.g. spin ladders and carbon nanotubes 
0, has been an additional source of motivation to investigate one- dimensional exactly solvable 
models. 

The basic concept of quantum integrability is the ^-matrix which represents either the 
factorized scattering of particles of (1+1) quantum field theories or the statistical weights of 
integrable two-dimensional lattice models. It turns out that the symmetry of the S'-matrix plays 
a fundamental role in the theory and classification of integrable systems H |J. Of particular 
interest is the Sp(2N) symmetry which preserves bilinear antisymmetric metrics, typical of 
systems with N component Dirac fermions. Even though the Bethe ansatz solution of the 
integrable Sp(2N) models has long been known |5|, ||, basic properties such as the nature of 
the elementary excitations and the low-energy behaviour have not yet been determined. 

The purpose of this paper is to unveil the physical content of the fundamental (vector 
representation) Sp(2N) solvable magnet. We argue that the low-energy properties are given 
in terms of 2N Majorana fermions due to the presence of special low-lying excitations in the 
spectrum. There exists at least two immediate applications of this result. First, it provides 
us a unique counterexample to the conjecture that integrable models based on the vector 
representation of Lie algebras should be lattice realizations of Wess-Zumino-Witten (WZW) 
conformal theories [^, || [|. In fact, we predict c = N for the fundamental integrable Sp(2N) 
model while the central charge of the Sp{2N) WZW theory is c = iV(2iV + 1)/(N + 2). Next, 
our study is of utility to one- dimensional systems with coupled spin and orbital degrees of 
freedom such as the spin-orbital [0, 11, [12|] and spin-tube |13| models. More precisely, we 
recall that the effective spin-isospin Hamiltonian describing these systems may be written in 
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the form |IT0|, |14 



H,so(Jo, Jl, J 2) — E E Ja^iM-l 
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where J Q are superexchange constants and i^ j+i denote the respective projections on the sin- 



glet, triplet and doublet spin-isospin states, see ref.jllj for details. By writing these projectors 
in terms of two commuting sets of Pauli matrices it is not difficult to identify that the in- 
tegrable Sp(4) spin chain [Hj corresponds to the point Jq/J\ = J0/J2 = 1/3. This point is 
interesting because it corresponds to both anisotropic and asymmetric spin- 1/2 couplings JTBJ ) 
being closer to represent the properties of realistic materials [19|, 2Q] than the integrable SU{4) 



case J = Ji = J 2 jig] . This then provides us a rare opportunity to determine exactly the 
nature of the excitations in a relevant spin-orbital model. 

In the context of statistical mechanics the integrable Sp(2N) model is a multistate vertex 
system defined on the square lattice whose bonds variables take 2N possible values. For each 
type of configuration of four bonds a, b, c, d meeting at a vertex we associate a Boltzmann 
weight factor S^(A) where A is the spectral parameter. Compatibility between integrability 
and the Sp(2N) invariance ("ice-type" restriction) leads us to the following amplitudes ||, || 
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A + iV + 1 

where a = 2N + 1 - a, e a = 1 for 1 < a < N and e a = -1 for iV + 1 < a < 2N. 

For any integrable vertex model one can associate a local spin chain commuting with the 
corresponding transfer-matrix whose matrix elements are given by ordered product of L factors 
S^(A). As usual, the Hamiltonian is proportional to the logarithmic derivative of the transfer 
matrix at the regular point A = 0, and in this case the expression is 
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where are the elementary matrix [e 



ab\l,k 



8 a ,i8b,k acting on site i. We observe that the 



spectrum of Hso{J\/^i Ji, Ji) matches that of JiHs p {a) ~ JiL. The Sp(2N) Hamiltonian (3) is 
solvable by the Bethe ansatz [|5], U and its eigenvalues E(L) can be parametrized in terms of a 
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set of variables A 



(a) 
j ' 



m a and a 



_J_ Va_ 

J ija 



■ S a , 1 



1, • • • , N, satisfying the following Bethe equations 
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and the eigenvalues are given by 
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where C a b is the Cartan matrix and r] a is the normalized length of the a-th root of the Sp(2N) 
algebra. 

We start our study considering first the Sp(4) model motivated by its direct relevance to 
the physics of spin-orbital systems. In fact, this is the simplest symplectic invariant system 
since N = 1 is equivalent to the isotropic six-vertex model. Later on we will show that the 
technicalities entering in the analysis of the Sp(4) model can be easily generalized to include 
arbitrary N > 2. Essential to our study is to determine the configurations of roots that describe 
the absolute ground state and the elementary excitations. This can be done by solving the 
Bethe equations (4,5) for some values of L and comparing it with the exact diagonalization of 
the Sp(A) spin chain. Its spectrum is parametrized by two set of variables A] ' and we found 
that the ground state and the low-lying excitations are characterized by strings with different 
lengths, namely 
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with 7 positive and C,^ a = 1,2,3 are real numbers. 

The important point here is to observe that the rapidities A^ 1 are described in terms of two 

independent type of strings, i.e 1-strings and 2-strings. This feature should be contrasted to the 

behaviour of the Bethe ansatz roots of other fundamental integrable systems such as SU(N) 

,(«) 



and O(N) models [|T6| , [8]]. In fact, for the latter systems each a-th root A!-" 7 characterizing 
the infinite volume properties has a unique string length. To explore the consequences of the 
peculiar string configuration (6) we substitute it in the Bethe ansatz (4) and we obtain the 
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following effective equations for the variables 
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where ip a (x) = 2 arctan(x/a), N a is the number of roots £j Q , Qy' = —(N a — l)/2+j — 1 with 
j = 1, • • • , iVo, and the matrix elements are given by 

<j>ap{x)= -lpl(x) 1p 2 (x) ~1pl/2(x) ~ ^3/2(x) (8) 

V tpl/2 (x) +*p3/2(x) -tpl/2 (x) ~ 1p3/2(x) 2^1 (x) + ^2 (x) J 

The ground state for even L consists of a sea of 1-strings and 2-strings with N± = N 2 = 



2N 3 = L/2. For large L, the roots £- a ^ are densely packed into its density distribution cr(^) a> ) 



1/L(£j"\ — ^ a> ) and the relations (7) in the L — > oo limit become integral equations for such 
densities. These integral equations are solved by elementary Fourier techniques and we find 
that 
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and 
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^ (3) 0r) = ^/ (rW(y)(7( 2 )(x-y)dy 
^= sinh(27ra;/3) — x 



6 sinh(7r:r) ^ 
where we have emphasized the remarkable fact that a^(x) is exactly the convolution of the 
densities a^\x) and a^ 2 \x). Recall that function a^{x) is related to the continuous proba- 
bility densities of finding the rapidity with a given value x. We may therefore interpret 
a^\x) as the probability for the sum of two independent events with probability cr^(y) at an 
arbitrary value y and a^(x — y) at the complementary value x — y. 

We now have the basic ingredients to investigate the thermodynamic limit properties. The 
ground state energy per site eoo is calculated by using Eqs. (5,9,10) after replacing the sum by 
an integral. The final result is 

I 

V3 ~ 3 
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The low-lying excitations are obtained by inserting holes in the density distribution of £j a ^ 
which means the removal of certain quantum numbers of the Bethe equations (7). The 
necessary manipulations of these equations are standard |T5| , [L6j and we find that the energy 



e ^(£) an d the momentum of one hole excitation in the sea of C,j a \ measured from the 

ground state, have the form 



eW(0=™ (a) (0» P (a) (0=y ? e {a \x)dx (12) 

where £ is the a-hole rapidity. For the first two excitations one can easily eliminate the variable 
£ leading us to the following dispersion relations 

e«(p) = - sin(2p), e< a >(p) = \ sin(2p) (13) 
2 o 

implying that these excitations are gapless and that their low-energy limits £^ a '(j>) ~ v ( a 'p are 
governed by distinct sound velocities, i.e = ir and = ir/3. The contribution to the 
total spin of each of these excitations is |, and therefore we shall interpret them as spinons 
propagating with different velocities. 

Similar computation for the third excitation leads us to transcendental equations and an 
analytical expression for the dispersion relation is hard to be obtained. However, it is possible 
to study the low-energy behaviour of such excitation by expanding the density (A) in powers 
of e~ A . We see that low-momenta regime is dominated by both sound velocities v^ 1 ' and 
and strictly in the p — > limit the lowest one prevails. This massless excitation turns out to 
be a spinless mode whose speed of sound is = n/3. At this point we note that recently the 
compound NaV^O^ has been modeled by an anisotropic/asymmetric spin-orbital model |I5 



Remarkably, the three-particle continuum found above is in accordance with the excitation 
spectrum proposed in ref.[|19[ to explain the optical properties of this material. 



Next we would like to identify the underlying conformal field theory which describes the 
low-energy limit of the integrable Sp(4) model. This can be investigated by analyzing the 
behaviour of the finite size spectrum of the corresponding spin chain. For a conformally 
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invariant theory, the ground state energy Eq(L) in a finite lattice of size L behaves as |2T|, |22 

!m =e .-£±w (14) 

where is the central charge associated to the a-th excitation. Similarly, the conformal 
dimension X- of the operator corresponding to the excited state Ei(L) is proportional to the 
finite size gap 



Ei(L) E (L) 2vr 
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Within the string hypothesis (6), the finite size corrections to the eigenspectrum can be 



evaluated analytically by applying the root density method [23] to the "string" Bethe equations 
(7). For instance, we find that each a-th massless excitation is associated with a central charge 

= 1. However, the string assumptions gives us only an idea of the true finite size behaviour, 
because the complex part of the roots may contribute to the term 1/L 2 as well. In order 
to obtain the correct finite size properties, we solve numerically the original Bethe ansatz 
equations (4) up to L = 36. In table (1) we exhibit such numerical results for the effective 
central charge c e f = Yh=i ^j c ^ an d the conformal dimension X e f = J2^=i associated 
with the lowest spin-wave excitation over the ground state. Clearly, the imaginary parts of 
the roots have conspired together and canceled the 1/L 2 correction proportional to the 
third mode. This is the first example we know that complex strings can produce negative 
contribution to the central charge. It is therefore tempting to think the third excitation as a 
composite state of two elementary spinons and that it does not contribute to the low-energy 
limit. The same phenomenon happens to the excited states and this allows us to interpret the 
operator content of the lowest excitation as X^ = x[ 2 ^ = 1/4. Putting these information 
together, we see that the underlying conformal field theory is likely to be represented in terms 
of four Majorana fermions rather than given by a Sp(A) WZW theory. 

Let us now turn to the problem of extending our results for general N > 2. The Bethe 
ansatz equations are parametrized by N different types of roots, and it turns out that the 
first N — 1 roots are given by both 1-strings and 2-strings while the last one behaves like 
1-strings. If we characterize the center of the strings by Q , a = 1, • • • , 2N — 1, we find that 
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the corresponding density distributions are 

( A. sm(na/N) _ _ _ AT _ 1 

p-W/^ — J ^ cosh(27rx/A0-cos(7m/A0 ' " ±; ' iV /jgN 

[ 2(AT+1) cosh[7rx/(Af+l)] ' a = N 

while for a = N+l, • ■ ■ , 2AT— 1 they are given by the convolution a^(x) = f^™ er( 27V ~ a ) (^o-W (x- 

y)dy. 

The corresponding excitations are gapless, consisting of N generalized spinons whose speed 
of sound are v^ a ' = ^ for a = 1, ■ • • , N — 1 and v = j^j- In addition, we have N — 1 
composite modes made by the convolution between the first N — 1 spinons with the iV-th 
excitation. For N > 2 numerical results for large L become difficult to be obtained since the 
number of roots to be determined grow rapidly with both iV and L. However, for N = 3 and 
small L ~ 18 our numerical analysis is consistent with the fact the only modes contributing 
to the low-energy properties are the spinons, each one with c = 1. All these results seem to 
be strong evidence that the continuum limit of such Sp(2N) integrable models can indeed be 
described in terms of 2N Majorana fermions. 

In conclusion, we have studied the excitation spectrum of the simplest integrable Sp(2N) 
spin chain. Contrary to the common belief this system is not the lattice realization of the 
Sp(2N) WZW conformal theory. Our study indicates that the nature of the excitations in 
spin-orbital systems can be rather involving. In fact, the isotropic point Jq = J\ — J2 is known 



to have three basic excitations [ 16| , being the lattice realization of a SU(4) WZW field theory 
0. However, the anisotropic point Jq/J\ = J0/J2 = 1/3 has only two independent excitations 
and one composite mode that do not contribute to the low-energy limit, and it is described by 
a c = 2 conformal field theory. This work prompts us to ask some questions that may open up 
new interesting avenues. What is the nature of the excitations of the spin-orbital model (1) in 
the crossover regime 1/3 < J0/J1 = Jq/ J2 < 1? What is the mechanism that made one of the 
excitations to become a composite state? What is the integrable lattice Sp(2N) model whose 
continuum limit corresponds to the Sp(2N) WZW theory? 



7 



Acknowledgements 

We thank A.L. Malvezzi for valuable discussions and the hospitality of the Institute of Ad- 
vanced Study, Princeton, where part of this work has been carried out. This work was sup- 
ported by Brazilian agencies Fapesp and CNPq. 

References 

[1] H.A. Bethe, Z.Phys. 71 (1931) 205 

[2] E. Dagotto and T.M. Rice, Science 271 (1996) 618; S. Iijima, Nature 354 (1991) 56 

[3] A.B. Zamolodchikov and ALB. Zamolodchikov, Ann. of Phys. 120 (1979) 253 

[4] M. Berg et al, Nucl.Phys.B 134 (1978) 125; V.V. Bazhanov, Phys.Lett.B 159 (1985) 321; 
M. Jimbo, Commun.Math.Phys. 102 (1986) 247 

[5] E. Ogievetsky and P. Wiegmann, Phys.Lett.B 168 (1896) 360; N.Yu. Reshetikhin, 
Lett. Math. Phys. 14 (1987) 125 

[6] For an algebraic Bethe ansatz derivation see, M.J. Martins and PB. Ramos, Nucl.Phys.B 
500 (1997) 579 

[7] I. Affleck, Nucl.Phys.B 265 (1986) 409; I. Affleck and F.D.M Haldane, Phys. Rev. B 36 
(1987) 5291; S.V. Prokovskii and A.M. Tsvelick, Sou. Phys. JETP 66 (1987) 1275; 

[8] M.J. Martins, Phys. Rev. Lett. 65 (1990) 2091; J.Phys.A:Math.Gen. 24 (1991) L159 

[9] A. Kuniba, Nucl.Phys.B 389 (1993) 209 

[10] K.I. Kugel and D.I. Khomskii, Sov.Phys.Usp. 25 (1982) 231 

[11] S.K. Pati, R.R.P. Singh, D.I. Khomskii, Phys. Rev. Lett. 81 (1998) 5406; Y.Q. Li, M. Ma, 
D.N. Shi and F.C Zhang, ibid. 81 (1998) 3527; B. Frischmuth, F. Mila and M. Troyer, 

8 



ibid 82 (1999) 835 ; P. Azaria, A.O. Gogolin, P. Lecheminant and A. A. Nersesyan, ibid 
83 (1999) 624 ! A.A. Nersesyan and A.M. Tsvelick, ibid 78 (1997) 3939 

[12] Y. Yamashita, N. Shibata, K. Ueda, Phys.Rev.B 58 (1998) 9114; S. Qin and I. Affleck, 
Phys.Rev.B 61 (2000) 6747; P. Azaria, E. Boulat and P. Lecheminant, ibid. 61 (2000) 
12112 



[13 
[14 

[is; 

[16 
[IT 
[18 

[19 

[20 

[21 

[22 



E. Orignac, R. Citro and N. Andrei, Phys.Rev.B 61 (2000) 11533 
D.P Arovas and A. Auerbach, Phys.Rev.B 52 (1995) 10114 

J.Des Cloizeaux and J.J. Pearson, Phys.Rev. 128 (1962) 2131; L.A. Takhtajan and L.D. 
Faddeev, Phys.Lett.A 85 (1981) 375 

V.G. Uimin, JETP Lett. 12 (1970) 225; B. Sutherland, Phys.Rev.B 12 (1975) 3795 
M.T. Batchelor et al, J.Phys.A:Math.Gen. 33 (2000) L97 



In the notation of ref. [|13l , equations (2.5-2.7), the Sp(4) point corresponds to two strongly 



coupled XXZ chains with J l = Jf = J 2 /2 = Jf = 2/3, A = 8/3, A 1 = A 2 /2 = 1. 

M. Mostovoy, J.Knoester, and D.I. Khomskii, \cond-mat/0009464\ (2000); M. Mostovoy 
and D. Khomskii, Solid State Commun. 113 (2000) 159 

S.K. Pati and R.R.P. Singh, Phys.Rev.B 61 (2000) 5868; A.K. Kolezhuk, H-J. Mikeska 
and U. Schollwock, \cond-mat/00lM2l\ (2000). 



H.Blote, J.L. Cardy and M.P Nightingale, Phys.Rev. Lett. 56 (1996) 742; I. Affleck, ibid. 
56 (1996) 746 

F. Woynarovich, J.Phys.A:Math.Gen 22 (1989) 427; A.G. Izergin, V.E. Korepin, and 
N.Yu. Reshetikhin, ibid 22 (1989) 2615 



9 



[23] H.J. de Vega and F. Woynarovich, Nucl.Phys.B 251 (1995) 439; F. Woynarovich and 
H.P. Eckle, J.Phys.A-.Math.Gen. 20 (1987) L97; C.J. Hamer, Q.R.W. Quispel and M.T. 
Batchelor, ibid 20 (1987) 5677; 



10 



J 


C e f 




1 o 


a n'v a rx^ a 
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1.000 985 
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1.004 993 


20 


4.U33 092 


1.006 011 


O/l 


/i noc; one 


1 nn^ ni7 


28 


4.020 376 


1.005 692 


32 


4.017 024 


1.005 272 


36 


4.014 620 


1.004 845 


Extr. 


4.001 (±2) 


1.006 (±2) 



Table 1 : Finite size sequences and the extrapolations of the effective central charge and the 
lowest conformal dimension of the Sp(A) spin chain. 
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